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In close analogy to the Bloch-Feshbach formalism known 
from the theory of nuclear dynamics, I develop a mathemati- 
cal framework that allows one to understand renormalization 
in terms of purely algebraic operations (projections, dilata- 
tions) in Hilbert space. This algebraic approach is put to the 
test in the study of the low-energy dynamics of interacting 
quantum gases, and proves to be efficient in deriving such 
diverse results as the renormalization group equation for an 
interacting Bose gas, the /3 function of (f>* theory, the screen- 
ing of fermion-fermion interactions or the BCS instability. 

ll.lO.Gh, 05.30.Ch, 71.10.-w, 21.30.Fe 



I. INTRODUCTION 

The idea of eliminating irrelevant modes in order to 
focus on the dynamics of few selected degrees of free- 
dom has a long history. The Bloch-Feshbach formahsm, 
developed in the latpn^O's to describe selected features 
of nuclear dynamics,a0 employs projection operators in 
Hilbert space in order to determine the effective Hamil- 
tonian in some restricted model space, thereby discard- 
ing dynamical information that pertains to the irrelevant 
modes. The irrelevant modes no longer appear explicitly 
in the effective theory, yet their residual influence on the 
dynamics of the remaining modes is taken into account 
through adjustments of the effective interaction. 

In a similar spirit, Wilson's renormalization group, de- 
veloped in the. early 70 's to better understand critical 
phenomena,^ Ll is a mathematical tool that allows one 
to iteratively eliminate short-wavelength modes and thus 
to arrive at effective ( "renormalized" ) theories which de- 
scribe the dynamics on successively larger length scales. 
In the original context of second-order phase transitions 
the renormalization group mainly served to calculate crit- 
ical exponents, and to provide a satisfactory theoreti- 
cal explanation for their universality. More recently, it 
has been pointed out that the renormalization group ra- 
tionale also affords a natural framework for the physics 
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of interacting fermion systems, and that it helps to un- 
derstand such diverse phenomena as Landau's fermi liq- 
uid theory, jj^arge-density waves, BCS instability, or 
screening.LTtHl 

The two calculational tools for systematic mode elim- 
ination -on the one hand the Bloch-Feshbach projec- 
tion technique in Hilbert space, based on algebraic con- 
cepts such as linear subspaces and operators, and on the 
other hand Wilson's renormalization group, commonly 
formulated in terms of functional integrals and Feynman 
diagrams- both derive their power from the fact that they 
allow one to study accurately selected features of the 
dynamics (for instance its infrared limit) without ever 
having to solve the full underlying microscopic theory. 
Beyond this common root for their success, mounting 
evidence like the success of Anderson^ "poor man's scal- 
ing" approach to the Kondo problem,llll Seke's projection- 
method treatment of the nonrelativistic Lamb shift Jia-ire- 
cent studies of a simple quantum mechanical model,ll3 or 
a calculation of the 1-loop renormalization oi-d)^ theory 
with the help of Bloch-Feshbach techniques,tJ suggests 
that the two methods are very closely related, and that 
in some cases the latter (Wilson) can even be regarded 
as a special case of the former (Bloch-Feshbach) . 

It is the purpose of the present paper to propose a 
slightly modified version of the Bloch-Feshbach formal- 
ism as an efhcient tool to do renormalization (Sec. ||). In 
the literature it has been put forth that the practical use 
of such a projection operator approach was limited to 
simple, essentially one-dimensional problems, and that 
the study of more complicated physical systems like a 
Fermi liquid would instead require the use of functional 
integral methodsfl Contrary to this assertion, I wish to 
show how algebraic techniques in Hilbert space can be 
successfully employed to tackle a variety of non-trivial 
phenomena in more than one dimension. I shall focus on 
the application to interacting quantum gases (Sec. Ill), 
both Bose and Fermi, and illustrate the versatility of the 
algebraic approach in deriving four well-known results: 
the flow of the coupling constant for bosons with point 
interaction; closely related, the one-loop /? function of 
0* theory; the screening of fermion-fermion interactions; 
and the BCS instability (Sec. 0). 

Despite the obvious parallels between the Bloch- 
Feshbach formalism and Wilson's renormalization, both 
in spirit and -as I will show- in practical calculations, 
there remains however an important difference: beyond 
the elimination of short- wavelength modes, which is com- 
mon to both approaches, the renormalization program 
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generally requires the additional step of rescaling. How 
this additional step, too, can be implemented within the 
algebraic framework, and under which circumstances it 
is called for, will be discussed in a separate section (Sec. 
0) at the end of the paper. 



II. PROJECTION TECHNIQUE 

We consider a macroscopic system (e. g., an inter- 
acting quantum gas) whose possible microstates span a 
-typically large- Hilbert space (e. g., boson or fermion 
Fock space), and whose dynamics in this Hilbert space 
is governed by some Hamiltonian H . In equilibrium the 
state of the system is described by a canonical statistical 
operator 



p{b) = exp[-i//6 - In Z{b)] , 6 > 



(1) 



with partition function Z{b), where b may represent the 
temperature T, Boltzmann's constant fc^, their product 
/csT, or possibly some other quantity which we choose 
to keep explicitly as an external parameter; it will be 
specified later. All other parameters, such as masses, 
coupling constants, the chemical potential or -in case 
b ^ {r, fcsT}- the temperature, are absorbed into the 
definition of the Hamiltonian. Given the statistical oper- 
ator as a function of 6, the underlying Hamiltonian can 
be extracted via 



b^^^p{b) = [H + c{b)]p{b) 



(2) 



Here c{b) is a c-number which stems from the deriva- 
tive of the partition function. Its specific value need not 
concern us, as long as we disregard simple shifts of the 
ground state energy. The partial derivative with respect 
to the external parameter is taken at fixed masses, cou- 
plings, chemical potential, and other parameters in the 
Hamiltonian. 

We presume that we are only interested in certain se- 
lected features of the macroscopic system (e. g., its long- 
wavelength properties), and that these selected features 
are represented by observables which merely act in some 
subspace of the original Hilbert space (e. g., in the sub- 
space spanned by many-particle states with all momenta 
below some cutoff). Let the operator which projects the 
original Hilbert space onto this selected subspace be de- 
noted by P, and its complement by Q = 1 — P. Informa- 
tion pertaining to the selected degrees of freedom iSj-then 
entirely encoded in the reduced statistical operatoitj 



p,s{b):=Pp{b)P ; 



(3) 



the other parts of the original density matrix, PpQ, 
QpP andyOpQ, only carry information which to us is 
irrelevant .Il3 

We would now like to cast the reduced statistical oper- 
ator again into the canonical form (|l|), with some modi- 
fied, effective Hamiltonian iJeff which acts in the smaller 



subspace only (ffcff = PH^sP)- This effective Hamilto- 
nian will in general acquire a dependence on the param- 
eter 6, so we must write 



Peff(fo) = P eM-H,s{b)/b ~ lnZeff(6)] 



(4) 



We assume, however, that the effective Hamiltonian is 
a smooth function of b which can be Taylor expanded 
around 6 = 0. For small parameter values we may then 
set, as a first approximation,E3 



i?eff(6) ~ i?cff(0) 



(5) 



This approximate effective Hamiltonian can be deter- 
mined in analogy to Eq. (||) , via 



= [i7eff(0)+C>eff(0) 



(6) 



&=0 



again up to some (generally different) c-number. 

In order to evaluate the left-hand side of the equation 
(^, we first use (||) to obtain 

6'^Pcff(6) - [PHP + c]peff (&) + PHQp{b)P (7) 

which includes a new term PHQp{b)P that accounts for 
the overlap between the P- and Q-sectors. Assuming 
that the original Hamiltonian can be decomposed into a 
free part and an interaction part, 



H = H^°'> + V 



(8) 



where the free part H^^^ commutes with the projection, 
this overlap term is 0{V^) and hence expected to be 
small for a weak interaction. With the help of the identity 



Qp{b)P = -Q j^J' dX ^ [e-^"Qe^"] p{b)P 



= -Q I dXe-^''[Q,V]e^"p{b)P (9) 
Jo 

we find, to lowest nontrivial order in perturbation theory, 
"exp(-£(")/6)-l 



Qpib)P = Q 



£(0) 



-iQVP) 



Pesib) . (10) 



Here denotes the Liouville "super" operator which, 
when acting on an arbitrary Hilbert space operator A, 
takes the commutator with the free Hamiltonian: 



(11) 



Since 



£(0) !£;(")) (ijW I ^ iEl"^-Ef)\El"^){Ef\ (12) 

for free energy eigenstates {\eI^^)}, it may be viewed as 
probing the energy difference between in- and out-states. 
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Provided the states in the Q-sector have higher ener- 
gies than those in the P-sector, the equation (|l^) has 
a well-defined limit 6^0: simply, exp(— £(°)/6) 0. 
Inserting this limit into Eq. (|^) then yields 

H^ff{Q) = PHP + J: + c" , (13) 

with S given to second order perturbation theory by 

= -pygJ^Q^p . (14) 

Clearly, the effective Hamiltonian is not just the projec- 
tion PHP of the original Hamiltonian, but contains also 
the extra term S to account for the residual influence of 
the eliminated Q-modes. 

III. APPLICATION TO INTERACTING 
QUANTUM GASES 

A. Mode Elimination 

We are interested in the low-temperature properties of 
interacting quantum gases, and hence in the effective dy- 
namics of low-energy excitations above the many-particle 
ground state. For noninteracting bosons the ground state 
has all particles in the lowest energy, zero momentum 
single-particle mode; while for noninteracting fermions 
the ground state consists of a filled Fermi sea, with all 
momentum modes occupied up to some Fermi momen- 
tum Kp. (For simplicity, the Fermi surface will always 
be taken to be spherical.) We will assume that, at least 
to a good approximation, these essential features of the 
ground state survive even in the presence of interaction. 
More specifically, we will assume that in the case of in- 
teracting bosons the ground state still has most particles 
in modes with zero, or at least very small, momentum; 
and that in the case of interacting fermions there still ex- 
ists a well-defined Fermi surface. Low-energy excitations 
then correspond to the promotion of bosons from small 
to some slightly higher momentum, or of fernuans from 
just below the Fermi surface to just above itxB At low 
temperature the regions of interest in momentum space 
are therefore the vicinity of the origin (bosons) or the 
vicinity of the Fermi surface (fermions), respectively. 

We now wish to devise a systematic procedure for fo- 
cusing onto these regions of interest. To this end we con- 
sider effective theories (i) in the bosonic case for modes 
within a sphere around the origin, of radius A; and (ii) 
in the fermionic case for modes within a shell inclosing 
the Fermi surface, of mean radius Kp and thickness 2A 
(where A <C Kp). Whereas in the limit of large A one re- 
covers the original, full theory, the opposite limit A — > 
yields the desired low-energy effective theory. In order 
to interpolate between these two limits we proceed in in- 
finitesimal steps: we lower the cutoff from some given 
A(s) to 



A(s + As) := exp(-As)A(s) , As > (15) 

with As infinitesimal, thereby discarding from the the- 
ory momentum modes pertaining to an infinitesimal shell 
(in the fermionic case: two shells) of thickness AA = 
A(s) As; determine the effective dynamics of the remain- 
ing modes; then eliminate the next infinitesimal shell, 
again determine the effective dynamics of the remaining 
modes, and so on. After each infinitesimal step we obtain 
a new effective Hamiltonian with slightly modified cou- 
pling constants. These may also include novel couplings 
which had not been present in the original theory; in 
fact, the mode elimination procedure will typically gen- 
erate an infinite number of such novel couplings. But 
in many cases only a few coupling constants will change 
appreciably and thus suffice to study the physical sys- 
tem at hand. How these coupling constants evolve as 
the flow parameter s increases and hence the cutoff A(s) 
approaches zero, can then be described by a small set of 
coupled differential equations. Modulo rescaling, which 
I shall discuss separately in Sec. these are the renor- 
malization group equations of the theory. 

Whether or not the temperature is among the quanti- 
ties that flow as a function of s, depends on whether it has 
been included in the definition of the external parameter 
h. If h is chosen to be ks or some other temperature- 
independent quantity then the temperature may flow, 
like all variables which have been absorbed into the def- 
inition of the Hamiltonian. If 6 = T or 6 = fc^T, on 
the other hand, then the temperature is regarded as a 
parameter that is controlled externally and hence fixed: 
it does not change upon mode elimination. While the 
former scenario applies to the study of isolated systems 
which, depending on the modes selected, may exhibit ef- 
fective dynamics at varying apparent temperatures, the 
latter scenario is adapted to the study of systems which 
are coupled to a heat bath of prescribed temperature. It 
is this latter case which we shall consider, as we wish 
to study the effective dynamics in the low-temperature 
limit T — > and thus explicitly control the temperature. 
Consequently, we choose h = ksT. 

Our aim is now to illustrate in a few prominent cases 
-bosons with point interaction, (j)^ theory, screening of 
fermion-fermion interactions, BCS instability- how the 
applicable renormalization group equations can be de- 
rived efficiently with the help of our projection tech- 
nique. To this end we must specify the appropriate 
many-particle Hilbert space (Fock space) and the appro- 
priate projection operator for each infinitesimal elimina- 
tion step. 

At a given cutoff A the boson Fock space is spanned 
by the particle- free vacuum |0b) and all n-particle states 
[n = 1 . . . oo) 

n 

|ki...k„) (x[]at(k,^)|0^) , |k,|<A , (16) 

i=l 

where the {k^} denote the particle momenta and {a^(ki)} 
the associated bosonic creation operators. 
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The fermion Fock space, on the other hand, is spanned 
by the fiUed Fermi sea (fermionic vacuum) |0t) and all its 
excitations which have particles above the Fermi surface 
and/or missing particles ("holes") below it, all within 
a shell of thickness 2A. In order to cast this into a 
mathematical formulation it is convenient to change co- 
ordinates, from the true particle momenta {K;} to little 
( "quasiparticle" ) momenta 



|K,;| -/<p)K,, 



and additional discrete labels 



(Ji := sign(|Kj| - Kp) 



(17) 



(18) 



This coordinate transformation K —t (k, a) is invertible 
except for modes which lie exactly on the Fermi surface. 
States above the Fermi surface are labeled a = 1, while 
those below are labeled cr = — 1. The allowed excitations 
in fermion Fock space then have the form (n = 1 . . .oo) 

n 

|k± . . . k±) cx l[[e{a,)a\k,,a,) + 0(-a,)a(-k„ a,)]|Of) , 
1=1 

|kf|<A , (19) 

where the {k^} denote the momenta of particles (+) 
or holes (— ), respectively, and {a^}, {a} the associated 
fermionic creation and annihilation operators. For sim- 
plicity, we have omitted any spin quantum numbers. 

It is now obvious which form the projection operator 
will have that is associated with the infinitesimal cutoff 
reduction (^5): if applied to any of the excitations ( p^ ) 
or (^|) it will simply multiply the respective state by a 

product of 9 functions, Yli ^(^ 
the new cutoff constraint. 



z^'\]4'^'>\), to enforce 



B. Hamiltonian 

We consider interacting quantum gases whose dynam- 
ics in the original, full Hilbert space is governed by a 
Hamiltonian of the form 



J J. 



ijkl 



(20) 



with kinetic energy H^*^^ and a two-body interaction V. 
By definition the single-particle energies include the 
chemical potential. The annihilation and creation oper- 
ators obey 



[ai, a]]zp = Si- 



(21) 



for bosons (upper sign) or fermions (lower sign), respec- 
tively. Each term in the Hamiltonian is normal ordered 
(:...:) with respect to the many-particle ground state. 



In the bosonic case this just coincides with the usual 
normal ordering: all annihilation operators to the right, 
all creation operators to the left. The explicit normal 
ordering of the Hamiltonian is then redundant. In the 
fermionic case, on the other hand, normal ordering means 
shuffling all operators which annihilate the fermionic vac- 
uum (ui with CTi = 1 or with aj ~ —1) to the right, 

all others {aj, with cr^, = 1 or a; with cr; = —1) to the 
left, thereby changing sign depending on the degree of 
the permutation. The explicit normal ordering of the 
Hamiltonian then becomes nontrivial. In both cases the 
normal ordering ensures that the energy of the (bosonic 
or fermionic) vacuum is set to zero, 

(Ob,f|ff|OKf) =0 . (22) 



C. Modification of the 2-Body Interaction 

As we discussed earlier, each mode elimination will 
yield an effective Hamiltonian that will generally contain 
a slightly altered mass, chemical potential, two-body in- 
teraction, etc., and possibly new interactions such as an 
effective three-body interaction. Here we shall restrict 
our attention to the modification of the two-body inter- 
action. In order to determine this modification we must 
consider the general formula ( p^ ) for the effective Hamil- 
tonian. As we have seen, the effective Hamiltonian is the 
sum of the projected original Hamiltonian, PHP, and an 
extra term S that accounts for the residual influence of 
the eliminated modes. For our purposes the projection 
PHP of the original Hamiltonian need not concern us, 
as it will not lead to any modification of the two-body in- 
teraction. Rather, we must investigate the ramifications 
of the extra term E. 

Application of the perturbative result ( p^ ) yields: 

abed ijkl 

xP -.ajalajai-. Q-^Q -.alalabaa-. P . (23) 

The two projectors P at both ends of the operator prod- 
uct ensure that all external momenta lie below the new, 
reduced cutoff; whereas the projectors Q in the center 
force at least one internal momentum to lie in the in- 
finitesimal shell which has just been eliminated. There- 
fore, at least one pair of field operators must pertain 
to the eliminated Q-modes, and hence be contracted: 
a(t)a(t) (Ob,f|a(t)a(t)|Ob,f). The product of the remain- 
ing six field operators can then be rearranged with the 
help of Wick's theorem, to yield a decomposition 



^ i](2) + j:]^> + s^^^ + Ej,^^ (24) 

(2) 

where each E„ is normal ordered and contains n field 
operators. The various terms shift the ground state en- 
ergy (n — 0), modify the mass, the chemical potential. 



.(2) 



.(2) 



.(2) 
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or more generally the form of the single-particle disper- 
sion relation (n = 2), modify the two-body interaction 
(n — 4), and generate a new effective three-body inter- 
action {n — 6). 

Since we want to focus on the modification of the two- 
body interaction we consider only the term with ri = 4. 
Calculating this term from Eq. (|2^) involves two con- 
tractions, at least one of which must pertain to the elim- 
inated Q-modes (see above), and both of which must go 
across the central Q{1/C^^^)Q (due to momentum con- 
servation). Neglecting the energy of the external modes, 
we find for bosons 



A{lk\V\ji)+ = -A 



and for fermions 



E or. \. {mV\ba) + {ab\V\jz)+ 

.ah 

(25) 



+ + (26) 



with 



^(ZS) 

lk\ji 



ah 



{<ya)9i-ab) - e{-aa)e{ab) 



its cross term 



x{la\V\bi)^{bk\V\ja) 



^(zs') ^ _^(ZS) 

lk\ji 



(27) 



(28) 



and 



(BCS) _ _ A rV^ 6{(Ta)9{ab) - 0{-(Ja)0{-(Tb) 

W - ^^Z. 2{ea + eb) 

ah 

x(lk\V\ba)^{ab\V\ji)-] . 



(29) 



The A in front of the sums signifies that at least one 
of the internal modes (a, b) must lie in the eliminated 
shell. In the bosonic case the modification of the two- 
body interaction can be associated with a 1-loop "ladder" 
diagram; in the fermionic case, on the other hand, there 
are three distinct contributions which, with hindsight, 
can be id£ntified with "zero sound" (ZS,ZS') and BCS 
diagramsJj The ZS contribution and its cross term ZS' 
account for particle- hole excitations (ctq = ±1, (Tb = =f1), 
while the BCS term describes 2-particle {da = (76 = +1) 
or 2-hole {aa = erf, = — 1) excitations. 



2U 

{lk\V\ji)+ = -^(5k.+k,.k,+k, 



(30) 



with the Kronecker symbol enforcing momentum con- 
servation, Q being the spatial volume, and U the cou- 
pling constant. Provided the magnitude of the external 
momenta k^, k^ is negligible compared to the cutoff A, 
momentum conservation implies that the internal modes 
a, b must both lie in the eliminated shell, and that hence 
ea — ^b — ^A- Application of the general formula ( |25| ) 
then yields 



AU = 



2rJeA 
where the sum 



A 



E <^k,+kj,k,+ki. 
|kaMkfc|<A 



(31) 



A 



E ^k.+kj,k„+ki, 
|ka|.|k6|<A 



A 



|ka|,|k6|<A 



E 

|k„|e[A-AA,A] 



(32) 



simply counts the number of eliminated states. For a 
spherical cut in momentum space this number of states 
is given by 



E 

|k„|e[A-AA,A] 



1 



(33) 



with p{e\) denoting the density of states at the cutoff. 
With AA = A As we thus obtain the fiow equation 



AU = - 



dine A p{€a ) 
rflnA 2n 



U'-As 



(34) 



For a dilute gas of nonrelativistic bosons in three spa- 
tial dimensions, with mass m, dispersion relation ea = 
A^/2m and density of states p(eA) = ri?nA/27r^ the flow 
equation reduces to 



AC/ = -^f/2.A. 
27r'^ 



(35) 



By its very definition the sequence of effective theories 
retains complete information about the system's low- 
energy dynamics. Observables pertaining to this low- 
energy dynamics are therefore unaffected by the succes- 
sive mode elimination, and hence independent of s. For 
example, the scattering lengthlij 



IV. EXAMPLES 

A. Bosons with point interaction 

For spinless bosons with point interaction {5 function 
potential in real space) it is 



m 
An 



U{s) - U[sf 



|p|<A(s) 



(Pp m 
(27r)3 ^ 



(36) 



stays constant under the flow (p5D, as the s-dependence 
of the parameters U and A just cancels out (up to third 
order corrections). 
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B. The link to 4>'*' theory 

There is an interesting relationship between the result 
(pj) and the 1-loop (3 function for real (j)^ theory. The 
(fr Hamiltonian describes the dynamics of coupled anhar- 
monic oscillators. It reads, in three spatial dimensions, 
H = + V with kinetic energy 



1 



(37) 



and interaction 



^ E n;^K. + «Uj^2..,- (38) 

4 Ct=l ^ 



kik2k3k 



Here m denotes the mass, n the spatial volume, g the 
coupling constant, and the single-particle energy 



ek 



(39) 



The field and its conjugate momentum tt are time- 
independent (Schrodingcr picture) operators which sat- 
isfy the commutation relations for bosons, and a, 
are the associated annihilation and creation operators. 
While the kinetic part of the Hamiltonian is normal or- 
dered (: ... :), the interaction is not. 

When expressed in terms of annihilation and creation 
operators the Hamiltonian takes on a form which is very 
similar to that of the quantum gas Hamiltonian (20) 



More precisely, the 0"* Hamiltonian contains a Bose gas 
Hamiltonian with two-body interaction matrix element 



3k,+k,,kfc-|-k, 



(40) 



The derivation of a flow equation for g can now proceed 
in the same vein as that for U, again starting from Eq. 
(^). Now, however, apart from 2-^2 particle scattering 
the Hamiltonian with its novel interactions a^a^a^a, 
a^a'^a^a^ etc. also permits 2 — > 4 and 2^6 scattering. 
Therefore, in Eq. ( p5| ) the intermediate state may be not 
just \ab), but also \abik), \abil), \abjk), \o-bjl) or \abijkl). 
As long as the magnitude of the external momenta is 
negligible compared to the cutoff, it is in all six cases 
Ea = Cfc = fA and 



{lk\V\ 
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{lk\V\jt 



(41) 



Hence in order to account for the larger set of allowed 
intermediate states we merely have to introduce an extra 
factor 6, and obtain thus 



A5 



d\n eA 3p(£a) 
din A Sflel 



g'-As 



(42) 



For A m it is eA = A, p{eA) = ^^a^^t^^, and the flow 
equation reduces to 



(43) 



in agreement with th o wel l-known 1-loop result for the /3 
function of (f)"^ thcory.E30 



C. Screening of fermion-fermion interactions 

We consider nonrelativistic fermions in spatial dimen- 
sion d {d > 2) which interact through a two-body inter- 
action 



{lk\V\jt)^ = [Viq)6,,sA.s, - V{ci')6s,s.S,,s, 



xS- 



Ki+Kj,Kfc+K, 



(44) 



duly antisymmetrized to account for Fermi statistics, and 
with {sq.} denoting the spin quantum numbers and q, q' 
the respective momentum transfers 



q 

q': 



Kfc 



K, = K, 



Kfc 



We investigate scattering processes for which 
0< |q|,A« |q'|,|K, + K,|,A> . 



(45) 



(46) 



In this regime only the ZS contribution (|27|) can signif- 
icantly modify the two-body interaction; its cross term 
ZS' (Eq. (H)), as wefl as the BCS contribution (|2|), 
are suppressed by a factor A/Kp. This can be seen 
directly from the geometry of the Fermi surface. The 
three constraints on the intermediate state: (i) both 
Ka and Kb lie in the cutoff shell of thickness 2A; (ii) 
more stringently, one of them lie in the infinitesimal shell 
to be eliminated; and (iii) — = — q' (ZS') or 
Kq, + Kft = Ki + Kj (BCS), respectively — reduce the 
momentum space volume available to the internal mo- 
mentum Ka to 0(Ar^^2AAA). In contrast, for |q| ^ A 
the ZS contribution with its condition K^ — K^ — — q 
allows a momentum space volume of the order AT^^^AA. 

To evaluate the ZS contribution at some given momen- 
tum transfer q, we first define the angle "i? between — q 
and the internal momentum K^, 



cos ■& 



K„ 



|q||Ka| 



(47) 



change coordinates from original (K) to little (k) mo- 
menta, and write, up to corrections of order |q|/A'i?, 

Ea - £& = W_f(|Kq| - |Kb|) = VF{<Ja\'i^a\ " O'bl'kbl) 

= VF\q\z (48) 
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with vp denoting the Fermi velocity. Next we note 
that the term with 0{aa)9{—ab) and the term with 
S{—cra)0{<Jb) yield identical contributions; therefore it 
suffices to consider only the first term and then multi- 
ply it by two. Finally, assuming that in the interaction 
matrix element the cross term is negligible, 



and thus 



|V^(q')l « l^(q)l 



(49) 



the two matrix elements in Eq. ( p7|) can simply be re- 
placed by V{q)'^ modulo Kronecker symbols for spin and 
momentum conservation. By virtue of these Kronecker 
symbols one of the two summations over internal modes 
collapses trivially, leaving 



AV{q) = A 

VF 



Y,e{A-\K\)e{\K]-\ci\z + A) 

^e{c7a)e{\q\z-\ka 



|q|^ 



(50) 



Here the first two functions explicitly enforce the sharp 
cutoff constraint for both ka and k;, (|ka|, |k{,| < A), while 
the latter two 9 functions enforce CTq = 1 and ab — —1, 
respectively. Under these constraints it is always G 
[0,7r/2) and hence z > 0. 

The above equation can be immediately integrated 
from cutoff A 3> |q|, |ka| (symbolically, A — > oo) down 
to A <C |q|, |ka| (symbolically, A 0), to yield the total 
modification of the two-body interaction: 



1 



1 



Kff(q) 

arc VH./ 



-|q|z-f A)^?(a,)0(|q|z-|kj) 



|q|^ 





(51) 



At the lower bound (A — > 0) the various conditions im- 
posed by the functions cannot all be satisfied simultane- 
ously, and therefore the product of functions vanishes. 
At the upper bound (A — > oo), on the other hand, the 
cutoff constraints imposed by the first two functions 
are trivially satisfied and thus can be omitted. In this 
case the sum over a is evaluated by turning it into two 
integrals, one over a radial variable such as |ka| or ea, 
the other over the solid angle. At a given solid angle and 
hence given z, the fourth function restricts the radial 
integration to the range |ka| € [0, |q|2;] or, equivalently, 
Ca G [0, Wii^lqjz]. This energy interval in turn corresponds 
to a number [p(ei;')wF|q|z] of states, pitp) being the den- 
sity of states at the Fermi surface. (It includes a factor to 
account for the spin degeneracy.) The integration over 
the solid angle is constrained to a semi-sphere, due to 
€ [0, 7r/2), and hence reduced by a factor 1/2 as com- 
pared to a full-sphere integration. Altogether we obtain 



1 



1 



K,ff(q) Vbare(q) 



- 2PM-F|qk^ 



(52) 



Kff(q) = 



Vbare(q) 



+ P(eF) 



(53) 



This result describes iie well-known screening of fermion- 
fermion interactions .t3 



D. BCS instability 

Our last example pertains to fermions with an attrac- 
tive pairing interaction 



{ik\v\ji) 



K„-K.OK, 



(54) 

which is the simplest form of BCS theoryii Due to the 
pairing condition Kj — — K^, K/ = — it is impossible 
to satisfy in the ZS and ZS' terms the requirement that at 
least one of the internal modes be in the eliminated shell. 
Hence only the BCS term ( p9| ) can modify the coupling 
constant. In the BCS term there are contributions with 
0{aa)0{<Jb) and 0{—aa)0{—(Jb), respectively, which yield 
identical results; therefore it suffices to consider only the 
first contribution and then multiply it by two. The pair- 
ing condition implies Cq = e;, = ca, which for modes in the 
upper (cr = +1) eliminated shell is given by ea = vpA. 
The eliminated shell itself covers an infinitesimal energy 
interval of width [wi?AAs], which in turn corresponds to 
a number [p(ei?)ui?AAs] of states. Of the two summa- 
tions over internal modes one collapses trivially due to 
momentum and spin conservation, leaving 



AV 



2vfA 
2vfA 



a 

piep) vpAAs 



V'-As 



(55) 



From this flow equation for the BCS coupling V{s) we 
immediately conclude that as long as the initial coupling 
V{Q) is positive, V{s) diverges as s ^ oo. This indicates 
the occurence of binding ("Cooper pairs") at very low 
temperatures. Furthermore, we can again convince our- 
selves that the sequence of effective theories retains com- 
plete information about the system's low-energy dynam- 
ics: Jew-energy observables such as the zero-temperature 
gapB 



Aq = 2A(s) exp 



p{ep)V{s) 



(56) 



do not depend on the flow parameter s and are thus un- 
affected by the successive mode elimination. 
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V. RESCALING 

So far we have only considered the systematic ehmina- 
tion of modes. Yet there is a second operation, namely 
rescaling, which in the context of renormalization plays 
an equally important role. In this section we shall first 
discuss the rescaling operation by itself, and later its in- 
terplay with the elimination of modes. 

Again we would like to couch the rescaling procedure 
into the language of linear operators in Hilbert space. 
To this end we define a unitary dilatation operator D{s) 
which, if acting on a (box-normalized) boson excitation 
( |l6| ) or fermion excitation simply rescales all mo- 

menta by e*:ca 



^?(s)|ki±^..k(±)) :=|e^ki±\..e^k(±)) Vs>0. (57) 

This operation can be viewed in two different ways, anal- 
ogous to the active and passive interpretations of coor- 
dinate transformations: either as an enlargement of the 
physical momenta at fixed length unit (active interpre- 
tation), or as an enlargement of the length unit at fixed 
physical momenta (passive interpretation). The latter 
view might be pictured as the observer with a televi- 
sion camera "zooming out" away from the probe, thereby 
watching on the camera's monitor a larger section of 
the probe and excitations with apparently smaller wave- 
lengths. 

Let the state of a physical system be described by a 
canonical statistical operator 



p{b) cx exp[-H{n)/b] 



(58) 



where we have written explicitly the possible dependence 
of the Hamiltonian on a finite volume fl. After "zoom- 
ing out" the canonical distribution will appear distorted, 
namely as 



/^zoom 

(6,s) ■.= D{s)p{b)D^is) 



(59) 



We wish to cast this zoomed distribution again into a 
canonical form, 



1 

(6, s) cx exp[-ffzoom(e ''Ti, s)/b{s)] 



(60) 



with a modified external parameter b{s) and some mod- 
ified, "zoomed" Hamiltonian -ffzoom- Since the rescaling 
makes the system appear smaller on the observer's (imag- 
inary) monitor, the modified Hamiltonian must now de- 
pend on the rescaled volume e~'^'''Ct {d being the spatial 
dimension) rather than on O. The thus defined zoomed 
Hamiltonian is easily determined: 

H,oon.{n,.s)^^-^D{.s)H{e+^-'n)D\s) . (61) 

In general, the s-dependence of the external parameter 
b is not a priori known and must be fixed by additional 



constraints. One such constraint might be the require- 
ment that some given reference Hamiltonian be scale 
invariant, 7J-f,„ = iJ-f. 

Just like the effective Hamiltonian obtained via mode 
elimination, the zoomed Hamiltonian generally contains 
modified masses, chemical potential and coupling con- 
stants; however, it does not contain any qualitatively new 
interactions. How the various parameters evolve as the 
flow parameter s increases, is again governed by a set of 
differential equations, the so-called scaling laws. 

To illustrate these general ideas we re-consider our flrst 
example, bosons with point interaction. Making use of 
the identity 



i:>(s)a(^)(k)L>t(s) =a(t)(e''k) 



(62) 



for annihilation and creation operators, we find for the 
free part of the Hamiltonian 



and for the interaction 



^oom 



^.(e-k) 



at(k)a(k) (63) 



-d-s 



bis) 



v{n) 



(64) 



Taking b = k^T, assuming a relativistic single-particle 
energy 



5(k) = ^/k^ + 



(65) 



with mass m and chemical potential /i, and requiring that 
in the limit m = fj, = the free dynamics become scale 
invariant (i?iSom = -H^*-"-*), we obtain the scaling relations 

T{s) = e'T , m(s) = e'm , ^i{s) = e'p (66) 

which in turn imply 

U{s)^e-'^'^~^'^'U . (67) 

This last result immediately translates into a scaling law 
for the coupling constant g of theory. The coupling 
g differs from U essentially by a factor ^eiCjek^i, which 
just leads to an extra factor e^* in the scaling: 



(68) 



As expected, the number e := 4— (rf+ 1) in the exponent 
is the deviation of the space-time dimension from the 
upper critical dimension four .123 

The combination of rescaling and mode elimination 
constitutes a complete renormalization group transfor- 
mation. In the by now familiar spirit we define the renor- 
malized statistical operator 



Pren{b,s) := D{s) P (s) p{b)P{s)D^ (s) 



(69) 
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which is obtained by first ehminating high momentum 
modes in the shell [e^^A,A] and then rescaling all mo- 
menta so as to recover the original cutoff. In our in- 
tuitive picture this transformation corresponds to first 
coarse-graining the physical spatial resolution and then 
"zooming out," thereby rescaling the length unit such 
that the apparent spatial resolution, defined in rescaled 
rather than original length units, stays constant. (One 
might visualize the apparent spatial resolution as the 
fixed pixel size of the observer's monitor, which, as the 
camera zooms out, corresponds to ever worse physical 
resolutions.) As in the previous cases we wish to cast the 
renormalized statistical operator into the canonical form 

p,en(^ s) cx exp[-i/rc„(e-''-^i^, bis),s)/b{s)] (70) 

with modified external parameter b{s) and renormal- 
ized Hamiltonian iJrcn- How the various parameters in 
the renormalized Hamiltonian (masses, couplings, etc.) 
evolve under the combined action of mode elimination 
and rescaling, is governed by the full renormalization 
group equations. These can be obtained by adding the 
respective flows due to mode elimination and rescaling. 

The rescaling part of the full renormalization group 
equations is specifically adapted to the study of critical 
phenomena. Critical phenomena are distinguished by the 
lack of an intrinsic length scale: as the observer is "zoom- 
ing out," at fixed apparent resolution, he keeps seeing the 
long- wavelength dynamics governed by the same masses 
and coupling constants. In other words, the observer has 
no means of inferring from these parameters the scale at 
which the dynamics is taking place. Critical phenomena 
are therefore naturally associated with a fixed point -the 
"critical" fixed point- of the full renormalization group 
equations. 

If one studies effects other than critical phenomena, 
such as screening or the BCS instability, then there is 
generally no reason for a repeated rescaling. 



VI. CONCLUSION 

In this paper we have formulated the renormalization 
program exclusively in terms of algebraic operations in 
Hilbert space. This new formulation has given us the op- 
portunity to expose the close relationship between Wil- 
son's renormalization and the projection technique a la 
Bloch and Feshbach, and to consider some of the ideas 
underlying renormalization from a slightly different con- 
ceptual perspective. Aside from these formal develop- 
ments we have also explored the potential for practical 
uses of the algebraic approach, and found that in four ex- 
amples -bosons with point interaction, 0** theory, screen- 
ing of fermion-fermion interactions, and BCS instability- 
the flow equations for the respective couplings could be 
derived in an efficient and straightforward manner. 

The algebraic approach carries the potential for inter- 
esting future developments. First, it is easily generalized 



since projection operators permit the elimination not just 
of short-wavelength modes, but also of other kinds of 
information deemed irrelevant, such as high angular mo- 
menta, spin degrees of freedom, or entire particle species. 
Secondly, the algebraic approach builds a bridge between 
renormalization and macroscopic transport theories, as 
also the transition from micro- to macrodynamics can ba 
effected by means of very similar projection techniques. E£l 
It may thus provide a natural language for the study of 
issues such as the renormalization of maccescopic trans- 
port theories, or effective kinetic theories.E^ 
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